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Abstract 

Based on our model of quantum systems as emerging from the coupled dynamics between oscillat- 
ing "bouncers" and the space-filling zero-point field, a sub-quantum account of nonlocal correlations 
is given. This is explicitly done for the example of the "double two-slit" variant of two-particle 
interferometry. However, it is also shown that the entanglement in two-particle interferometry is 
only a natural consequence of the fact that already a "single" two-slit experiment can be described 
on a sub-quantum level with the aid of "entangling currents" of a generally nonlocal nature. 
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1. INTRODUCTION 



Although nonlocality has featured very prominently throughout the last decades in the 
discussions on the foundations of quantum mechanics, no general consensus has yet been 
reached over it. Despite many arguments in favor of the position that nonlocal effects are one 
of the (if not the) main characteristics which distinguish classical from quantum mechanics, 
some researchers even hold on to the view that a purely local physics could suffice to explain 
all existing experimental data, thus ultimately relying on a familiar trait of classical physics. 

However, even if one considered quantum mechanical nonlocality as a well-established 
fact, as we do, this does not necessarily mean that there cannot exist some form of "classical" 
explanation for it. In arguing for the use of modern, "21 st century classical physics", our 
group, for example, has in recent years obtained with such classical means a series of results 
that were previously considered as obtainable only via quantum mechanics. Among these 
results, some features figure prominently, like, e.g., explanations of Planck's postulate of 
energy quantization, the dispersion of a Gaussian, or interference at a double slit. (For an 
introduction, see [1] and the references therein.) Thus, as a further task along our lines of 
reasoning, we intend in this paper to provide a classical framework for nonlocal effects and 
entanglement. 

One of our main modeling scenarios is provided by the stimulating experiments of 
Couder's group [2-6], where "bouncing", and also "walking", droplets (the "particles") are 
dynamically coupled to the oscillations of a bath (the "waves") and thus produce a whole 
series of effects, which were previously considered only to be possible quantum mechanically. 
Among these effects, interference at a double slit, tunneling, or quantization of angular mo- 
mentum in closed orbits could be reproduced, for example. Although it is clear that the 
mentioned experiments can only provide analogies, at best, one has here, nevertheless, a 
scenario providing essential stimuli for model building also in the context of quantum theory. 
This, at least, is what we want to propose here, i.e., that there are further insights to be 
gained from the experiments of Couder's group, which could analogously be transferred into 
the modeling of quantum behavior. Concretely, we do believe that also an understanding 
of nonlocality and entanglement can profit from the study of said experiments. In fact, one 
indispensable prerequisite for these experiments to work, one basic commonality of all of 
them, is that the bath is vibrating itself. It so happens that bouncer and bath may engage 
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in a self-organized, dynamically coupled entity where the bouncer can self-propel due to its 
interaction with the wave it generates. 

In recent papers, we introduced an analogous scenario in the quantum domain, thereby 
having a particle/bouncer undergo also stochastic jumps such that diffusion theory can be 
applied. Then, what is in a quantum mechanical context described as a Gaussian wave- 
packet has an equivalent in the Gaussian distribution of a particle whose path follows the 
"agitations", or "excitations", of the underlying "bath". With the latter, we refer in our 
model to the zero-point oscillations of the vacuum, i.e., something we take as empirically 
given, or, in other words, as an "ontological" input to our theory. Note that we completely 
agree here with Timothy Boyer who noted in a similar context: 

"The concept of zero-point radiation (random radiation fluctuations at the zero of 
temperature) can appear in both classical and quantum theories. Zero-point radiation 
can not be regarded as belonging exclusively to quantum theory any more than the 
concepts of mass, energy, and gravity can be claimed as exclusively classical concepts 
because they appeared first in the context of classical mechanics." [7] 
As our model's particle through its bouncing and locking-in with the zero-point field 
creates diffusion wave fields, and as the latter, at least in the nonrelativistic case, extend 
instantaneously along nonlocal distances [8-10], we speak of a nonlocal "path excitation 
field". Wherever the diffusion wave fields "radiated out" into the environment exist (or 
better: Wherever the bouncers' oscillations lock in with the oscillations of the space-filling 
bath), the self-propelled bouncer (= "walker") may go, eventually. So, in the case of a 
bouncer oscillating in some source region of an interferometer experiment, whenever it is 
propelled forward and can potentially go through one of the two slits, the (Gaussian) path 
excitation fields behind the double-slit will overlap to produce the familiar fringes at some 
screen [11]. We stress, however, that one of the most important features of the path exci- 
tation field is its nonlocal nature. This is actually what we are going to work out in more 
detail in the present paper. 

Apart from accepting as given the nonrelativistic diffusion wave fields' "breathing" non- 
locally throughout space, what could be a possible origin of the assumed nonlocal nature of 
the zero-point oscillations? One can only speculate at present, but it seems that a good can- 
didate for an explanation would come from cosmological considerations. Note, for example, 
that for the universe in its initial phases, according to present-day models, one can admit, in 
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addition to the particles existing in the very early universe, a set of phase-locked wave-like 
oscillations that would thus "resonate" throughout the whole small-scale universe. Then, it 
is conceivable that cosmic inflation, for example, would not destroy these oscillations, but 
rather "inflate" these fields as well, thus ending up with a much larger universe where the 
particles still oscillate in phase with the zero-point background, albeit with the latter now 
having turned a nonlocal one. Whatever the true reason for the nonlocality of the zero-point 
field may eventually be, we take the latter as an input for our modeling, and in the following 
show some results thereof. 

2. A SUB-QUANTUM KINEMATIC ACCOUNT OF TWO-PARTICLE CORRE- 
LATIONS 

Let us first consider an EPR-type experiment as it was proposed by Home and Zeilinger [12] 
for two-particle interferometry, i.e., essentially a "doubled" double-slit experiment as shown 
in Fig. 2.1. One creates a particle pair with equal but opposite momenta ki and k2, respec- 
tively, along with some small changes 5k, which can be varied via phase shifters $i and $2- 
Their quantum mechanical wave function would be described in terms of an entangled state. 
The resulting correlated intensity I (x 1; x 2 ) exhibits marked modulations, which persist even 
for arbitrarily large spatial separations of the individual particles' wave packets. 




Figure 2.1. Scheme of a two-particle interferometer, with the source S in the center emitting two 
anti-correlated particles, with different phase shifters 3> inserted into the particle beams, and with 
(possible locations of) detectors D. 

In our approach, we can make use of the path excitation field as follows. Each path 
% be occupied by a Gaussian wave packet with a "forward" momentum pj = Kki = mv,. 
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(Moreover, due to the stochastic process of path excitation, the latter is represented also by 
a large number of consecutive Brownian shifts, p UjQ = mu a , but we shall return to this only 
in the next chapter, where we shall discuss in more detail the different roles played by the 
velocity fields v and u, respectively.) To start with, we note first the probability density 
currents for each "particle", and then combine them in a suitable manner. Thus, upon 
preparation of the two-particle source, we start with some initial distribution P (xo, to) — 
Rq (x , to) °f a composite system just at decay time t . Considered as a two-particle system, 
and since we employ a wave theory, this initial distribution can be rewritten as P (x , t ) = 
Rl(x , t ) Rr (x , t ), with R L and Rr signifying normalized amplitudes associated with 
particles going left (locations xi) or right (mirror locations X2), respectively. This initial 
distribution is thus split up into one channel for particle 1 and a correspondingly anti- 
correlated channel for particle 2 (Fig. 2.1). Note that, because each of the two bouncers 
"excites" the areas in both directions (i.e., both possible paths) of the surrounding medium, 
the probability density for each path is given by P (x i; t) = R L (x i; t) Rr (x i; t), for i — 1 or 
2, and similarly for the primed quantities, where Rl and Rr now more specifically refer to the 
amplitudes "excited" by the bouncer going left or right, respectively, irrespective of which 
actual particle one focuses upon. For simplicity, we shall in the examples below concentrate 
on the symmetric scenario with equal weights, i.e., Rl (xi) = Rr (xi) = Rl (x 2 ) = Rr (x 2 ). 

Then, in close similarity to our treatment of the case of a single double-slit [ ], for the 
"doubled" setup we write down the total average probability current as the sum of all four 
average probability currents present, i.e., firstly without the presence of any phase shifters 
it holds (with bars denoting averages) that 



'tot 
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(2.1) 



Allowing now for a relative phase 7? = $1 — $2 through the insertion of the phase shifters $1 
and $2, one can write with the "total" average momenta k := ki + k2 and k := k : + k 2 + <5k, 
respectively, 
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RlRr (k : + k 2 ) + R' L R' R (k[ + k' 2 + 5k) 



(2.2) 



Then we obtain with normalization A/", with the momentum balance k tot = k = k', and 
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with hats denoting average unit vectors, the correlated intensity 



/ ( Xl , x 2 ) := APP t 2 ot (xi, x 2 ) = Af R L R R k + R' L R' R k' 



(2.3) 



and thus 



/(x 1)X2 ) =M 2 



p2 p2 I p/2 p/2 



+ 2R L R R R' L R' R cos | (ki + k 2 ) - (ki + k' 2 + 5k) ■ r} 



(2.4) 



This is the exact quantum mechanical result, albeit here obtained without invoking the 
quantum mechanical calculus. Moreover, one can now also highlight particular features 
of this remarkable correlation along nonlocal distances r = x\ + x 2 between the locations 
X! and x 2 , respectively. Namely, as opposed to the formula (2.4) derived from the total 
probability density current, one can also single out correlations between individual currents, 
respectively, for the various channels. An important role is thereby played by the relative 
phase if, which emerges naturally in our model as it relates a "bouncer's" oscillations along 
different paths [11, 13]. 



As is usual in interferometry, differences in relative phase are - also classically - accounted 
for by phase shifts of | for each reflection of a beam at one of the slabs of the interferometer. 
Thus, by comparing primed versus unprimed scenarios, one can for example relate the 
intensities at detectors D 2 and -D4 in Fig. 2.1, which provide the conditional probability 



P{D 2 \D,)=M 2 R 2 L R 2 R 
1 



2 + 2 cos I (ki + k 2 ) - (ki + k' 2 + 5k) ■ rj 



4 (2 + 2 cos [$!-$ 2 ]) = -(1 + cosp) 



(2.5) 



whereas 



P (£> 2 I D 3 ) = - (2 + 2 cos [$! - $ 2 + tt]) = - (1 - cos ip) . (2.6) 

However, 
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P( J D 6 | J D 4 ) = i(2 + 2cos£) =\, (2.7) 

for example, because the phase difference between the two possible paths is now independent 
of the $j and given by if = |. From the latter example one sees that an "early" detection of 
a particle at D$, which is equivalent to a which- way measurement, or to the "closing of the 
second slit" on one side of a double-double-slit experiment, respectively, destroys nonlocal 
interference effects such as those indicated by (2.5) or (2.6). 

How can we understand this behavior in our model? Now, we have mentioned on several 
occasions that we consider the Gaussians employed in our calculations as simple solutions of a 
diffusion equation. More generally, however, as has been explicated in the works of Mandelis 
et al. [8-10], diffusion wave fields related to oscillating sources may in a more specific way 
extend nonlocally across the whole domain of an experimental setup, for example. Thus, 
the Gaussians used so far may be only approximations to more complex solutions of the 
diffusion equation, which actually may exhibit long wiggly tails, albeit with very small 
amplitudes. However, it is exactly such a functional characteristic which has been found 
also in a specific quantum mechanical context, i.e., for example, in Rauch's post-selection 
experiments [14, 15]. There, it turned out that for interference to occur it is not necessary 
that the "main bulks" (approximated by Gaussians) of wave-packets overlap. Rather, the 
experiments show that interference can be caused by the nonlocally far-reaching action of 
the plane- waves of a quantum mechanical wave-function. This is actually the corollary of our 
understanding of the nonlocal nature of the zero-point field to which our particle/oscillator 
couples: A Gaussian indicating the approximate whereabouts of our bouncer, embedded in 
an oscillatory "bath" of momentum fluctuations, the regular part of which thereby coinciding 
with the action of the plane-wave components in quantum mechanics. 

Note particularly that the bandwidth of these plane-wave components is determined by 
the momentum resolution of the whole measurement apparatus, where the upper limit is 
defined by the inverse of the distance between source and detector [15]. This is strongly 
reminiscent of our analogy with the particle's bounces in the Couder experiments, which 
are locked-in with the oscillations of the fluid. The spatial constraints of the latter, i.e., 
the container sizes, thus define the momentum resolution of the experiment via a suitable 
bandwidth of possible wavelengths. So, we observe that in our model it is the wholeness 
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of the possible wave configurations within the possibly nonlocal limits of an experimental 
setup that co-determines the experiments' outcomes. 

In other words, whenever the "constraints" of the experimental setup are changed, this 
may have a nonlocal effect on the registered particles - a possibility which brings us back 
to the effect of the "closing of a slit" in the double-slit experiment. Note that, as this effect 
is shown in a forthcoming paper [13] to be essentially nonlocal, it is of the same nature as 
the one in EPR-type experiments, as, for example, in Aspect's experiments: As opposed 
to mere kinematic nonlocality like the one implied by intensity correlations such as (2.4), 
we are in ref. [13] interested in the effects of the actual "closing of a slit", i.e., in dynamic 
nonlocality [16]. However, in the present paper we shall restrict our discussion to the simple 
double-slit experiment. For this alone already suffices to bring forth essential features of 
nonlocality and entanglement: The classical roots of entanglement are already visible in the 
"single" double-slit experiment, as will be shown now. 



3. ENTANGLING CURRENTS IN THE DOUBLE-SLIT EXPERIMENT 

We consider the particles as emerging from one of two "Gaussian slits", i.e., two possible 
paths of a particle which later eventually cross each other. To do so, we started in [11] 
with an independent numerical computation of two Gaussian wave packets, with the total 
distribution given by 



Ptot :=Pi + P 2 + 2 v / ^cos^ 12 , (3.1) 

where the phase difference 

''- / (x - X 02 - V 2 t) 2 - (x - Xqx - Vit) 2 ' 



1 

n 



m{v2 _ Vl)x ~ ^ ~ ^' ) t 



(3.2) 



is characterized by the usual "classical" velocity difference vi — wiand a kinetic energy term 
including a momentum fluctuation muo, or the "osmotic" velocity uq, respectively. Here, 
we indicate the two slits at positions Xqi and £02 and we assume the same slit widths and 
hence the same initial standard deviations ctq. The graphical result of a classical computer 
simulation of the interference pattern in a double-slit experiment, including the average 
trajectories, with evolution from bottom to top, is shown in Fig. 3.1. The Gaussian wave 
packets characterized by moderate spreading at the same standard deviations a move to- 
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Figure 3.1. Classical computer simulation of the interference pattern in a double-slit experiment; 
with v Xj i = -v x>2 . 




Figure 3.2. Same as Fig. 3.1, with an additional phase A<p = tt at slit 1. 



9 



t 




XQ1 Xq 2 ^ 

Figure 3.3. Classical computer simulation of the total average "entangling current" density in a 
double-slit experiment; same setup as in Fig. 3.1, with arbitrary normalization and v x \ = —v x i- 




Figure 3.4. Same as Fig. 3.3, with an additional phase shift of Aip = tt at slit 1 according to 
the setup of Fig. 3. 2. Comparing with Fig. 3.3, one notes that the dramatic shift from maxima to 
minima, and vice versa, as observed in the interference patterns of Fig. 3.1 and Fig. 3. 2, respectively, 
is essentially caused by the changes in these entangling currents. 
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wards each other with velocities v Xi \ = —v x ^- One can observe a basic characteristic of 
the averaged particle trajectories, which, only because of the averaging, are identical with 
the Bohmian trajectories. To fully appreciate this characteristic, we remind the reader of 
the severe criticism of Bohmian trajectories as put forward by Scully and others (see [17], 
and references therein). We can note that in our sub-quantum approach an explanation of 
the obvious "no crossing rule" is straightforward and actually a consequence of a detailed 
microscopic momentum conservation. In Fig. 3.1 the maximum of the resulting distribution 
is positioned along the central symmetry line in between the two slits. [1 1] 

The interference hyperbolas for the maxima characterize the regions where the phase dif- 
ference (f\2 = 2mr, and those with the minima lie at ipi 2 = (2n + 1)tt, n = 0, ±1, ±2, . . . Note 
in particular the "kinks" of trajectories moving from the center-oriented side of one relative 
maximum to cross over to join more central (relative) maxima. In our classical explanation 
of double slit interference, a detailed "micro-causal" account of the corresponding kinemat- 
ics is given. (For the details, see [11].) Since each Gaussian has its own phase distribution 
(3.2), we are free to add a phase shift Aip at one of the two slits, say slit 1, which modifies 
(fi to 

<p l = ^- + A(p. (3.3) 
n 

In Fig. 3.2, we use the same double-slit arrangement as in Fig. 3.1, but now include a 
phase shifter affecting slit 1, as sketched by the yellow rectangle on the left hand side. 
The exemplary choice of Aip = it results in a shift of the interference fringes. Comparing 
with Fig. 3.1, we recognize now a minimum of the resulting distribution along the central 
symmetry line. 

Finally, we reconsider our classically obtained total average probability current [1, 11] 

■Jtot = P\V\ + P 2 v 2 + yPiP~2 {vi + v 2 ) cos V?12 + \JP\P 2 {u\ - u 2 ) sin if 12 , (3.4) 

with osmotic velocities Ui and convective velocities applied to both slits, i = 1 and 2, and 
with the phases (3.2) and (3.3). 

The result of our computer simulation of Eq. (3.4) is shown in Figs. 3.3 and 3.4 corre- 
sponding to the intensity distributions of Figs. 3.1 and 3.2, respectively. One recognizes the 
change of the maximum values of the probability current along the central symmetry line in 
Fig. 3.4 in comparison with those of Fig. 3.3. Since the figures display the one-dimensional 
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case, the current flow is along the x-axis only. Interestingly, at the time t T of the reversal 
of the trajectories, the current flow comes to a hold, and starts again for times t > t T with 
reversed signs. This can be understood as a reversal of the relative flow of heat Q2 — Qi 
between the two channels, since Ui = -^VQi [18], such that the last term of Eq. (3.4) 
reads as -^y/PxKV {Q 2 - Qi) sin (p 12 . 

The probability current Jtot i n both figures essentially only consists of the last term of 
Eq. (3.4), as the velocities Vi and the velocities typically differ by many orders of mag- 
nitude. In other words, the probability current J tot is always dominated by the "quantum 
mechanical" entangling term of Eq. (3.4) which is connected to the osmotic velocities, u\ 
and 1*2, and which implies the existence of nonlocal correlations. As we have just seen, this 
"entangling current" can also be understood as describing the "heat flow" between the two 
channels. Note that, as opposed to the average total probability current J tot , in the distri- 
bution (3.1) of the probability density P tot alone, just as in the correlated intensity (2.4) of 
two-particle interferometry, the entangling part is not explicitly visible. 

The phenomenon of entanglement is thus possibly rooted in the existence of the path 
excitation field, a version of which we already encountered in the present paper with the 
usual double-slit interference. 

In other words, one can say that the entanglement characteristic for two-particle in- 
terferometry is a natural consequence of the fact demonstrated here, i.e., that already in 
single-particle interferometry one deals with entangling currents, which generally are of a 
nonlocal nature. 
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